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The photoproduction of vector mesons on various nuclei has been studied using the CEBAF 
Large Acceptance Spectrometer (CLAS) at Jefferson Laboratory. The vector mesons, p, uo, and 0, 
are observed via their decay to e'^e~ , in order to reduce the effects of final state interactions in 
the nucleus. Of particular interest are possible in-medium effects on the properties of the p meson. 
The p spectral function is extracted from the data on various nuclei, carbon, iron, and titanium, 
and compared to the spectrum from liquid deuterium, which is relatively free of nuclear effects. We 
observe no significant mass shift for the p meson; however, there is some widening of the resonance 
in titanium and iron, which is consistent with expected coUisional broadening. 
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The relation between the medium modification of the masses and chiral symmetry restoration in high density and 
temperature has been the subject of many studies, both experimentally and theoretically, for about two decades. 
Spontaneous breaking of the chiral symmetry leads to a non-zero value of chiral condensate, {qq) ~ 300 MeV, which 
is responsible for the masses of hadrons Indeed, most of the mass is generated dynamically. As the value of {qq) is 
predicted to decrease with increasing ternperature and density, in-medium modification of the masses are expected in 
hot/dense mediums 0, 0, 0, S & 0, 11, S 113, El • Quantitative understanding of the dynamics in the non-perturbative 
regime of QCD is rather incomplete and therefore experimental data is essential. If these modifications exist, they 
could also be an indication of a transition from hadronic matter to a deconfined quark-gluon plasma [12]. 

To investigate these modifications, high temperature and density can be achieved by experiments using heavy-ion 
colliders. The first indication of a possible medium modification of the p meson came from the CERES [13] and 
HELIOS/3 [13] collaborations at CERN in 1995. CERES reported on the measurement of low mass e+e^ pairs from 
p-Au and Pb-Au collisions. While their proton-induced data could satisfactorily be accounted for by summing various 
hadron decay contributions, an enhancement over the hadronic contributions was observed for the Pb-Au data in the 
mass range between 300 and 700 MeV. However, the resolution and statistics of the data did not allow a quantitative 
measure of a mass shift or a change in the width of the p vector meson. Recently, following an upgraded experiment, 
CERES reported a broadening in the mass of the p vector meson due to the baryonic interaction in the low mass 
region [l^. The result of di-muon measurement in In-In collision in NA60 experiment at CERN has also shown a 
considerable broadening (doubling of the width) of the p spectrum, while no shift in the mass was observed [l^, [13, [3 ■ 

Although a large medium effect is expected in the heavy-ion reactions, the results are integrated over a wide range 
of temperature and density, many channels contribute, and the reactions proceed far from equilibrium. This makes 
an exact interpretation of results difficult, and the connection between chiral symmetry restoration and in-medium 
modifications remains unclear [Tol. [20l, [2l| . However, the in-medium effects for the vector mesons at normal nuclear 
density and zero temperature are predicted to be so large that they can be studied by fixed-target experiments 
involving elementary reactions [2^ . 

Brown and Rho Q start from an effective Lagrangian at low energy and zero density. They apply the same 
Lagrangian at nuclear density but with masses and coupling constants that are modified according to the symmetry 
constraints of QCD (e.g., chiral symmetry). This model predicts an in-medium scaling law that results in a decrease 
in the mass of vector mesons by 20%. 

Hatsuda and Lee based on QCD sum rule calculations, obtain mass changes of the vector mesons in the nuclear 
medium. Their calculations result in a linear decrease of the masses as a function of density introducing a mass shift 
parameter a: 

— l~a — , a = 0.16 ± 0.06, (1) 



mvAiip^O) po' 

where mvM is the mass of the vector meson, pQ indicates nominal nuclear density (0.16 fm^'^), and p = indicates 
the vacuum. 

Models based on nuclear many-body effects predict a broadening in the width of the p meson with increasing 
density. This prediction is based on the assumption that many-body excitations may be present with the same 
quantum numbers and can be mixed with the hadronic states [l^, |Tl|, [l^, lU, [2^, [2^. Furthermore, due to the 
uncertainty of the coupling constants as a function of density, the branching ratios are expected to change in the 
nuclear medium and also distort the invariant mass spectrum of the resonance (27j . 

An observation of a medium-modified vector meson invariant mass decrease {a = 0.09±0.002) has been claimed 
by a KEK-PS collaboration in an experiment where 12 GeV protons were incident on nuclear targets, and the e+e^ 
pairs were detected '^^H^, 113]. Very recently, the Crystal Barrel/TAPS collaboration has reported a 14% downward 
shift in the mass of the w, where the analysis focused on the ir'^j decay of low-momentum uj mesons photoproduced 
on a nuclear target [Slj . 

The data for the present study were taken in 2002 using the CEBAF accelerator and the CLAS detector located in 
Hall-B of Jefferson Laboratory. A comprehensive description of CEBAF, the Continuous Electron Beam Accelerator 
Facility, can be found in the Ref. [13], while Ref. [11] gives a detailed account of CLAS, the CEBAF Large Acceptance 
Spectrometer, and other Hall-B equipment. The incident bremsstrahlung photon beam on the target was produced 
from a primary electron beam of 3 GeV in energy for the first 2/3 of the run time, and a primary electron beam of 4 
GeV for the last 1/3 of the run time. 

The Cerenkov counters (CC) in combination with the electromagnetic calorimeters (EC) were the two most crucial 
CLAS components for this experiment. The EC and CC cover the forward part of the CLAS detector, subtending 
scattering angles of 8° < < 45°. The e+e~ event selection and the rejection of the very large tt+tt" background 
were done through cuts on the EC and the CC. The pion rejection (or misidentification) factor is determined to be of 
the order of 10^^ for a two-track measurement 1341. 
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The target contained a liquid-deuterium (LD2) cell and six solid foils, each with a 1.2 cm diameter. Four of the foils 
were carbon, one was titanium, and one iron. The total thickness of the deuterium, and the four carbon targets was 
each 1 g/cm^, while the titanium and iron targets were each ^ g/cm^. The atomic weights of iron and titanium were 
close enough that the data from these two targets were combined to increase the statistics. The separation between 
targets was 2.5 cm; the target nucleus was determined by the reconstructed position of the production vertex; and 
the CLAS vertex reconstruction resolution was 0.3 cm. 

The object of this study is the invariant mass of e^e~ from the decay of the vector mesons p, w, and (p. This 
branching ratio for vector mesons into e+e~ is of the order of 10~^ — 10~^. Other physical processes also contribute to 
the background, for example lu ir^e^e^ , rj 76+6", and n'^ je^e" (Dalitz decay). In the case of the 77 and 7r°, 
the e~^e~ mass is below the region of interest. In the case of the lu, the Dalitz decay is also included in the expected 
spectrum, with this mode tied to the e~^e~ mode. The background from — > tt^tt^X with both pions decaying 
via the Dalitz mode is also considered. In this case, the e+ may be detected from one pion and the e~ detected 
from the other. This process was simulated with the known n^n^ cross section (35| and angular distributions, and 
its contribution determined to be small (0.02%). In addition, Bethe-Heitler e+e" pairs have been simulated with the 
expected cross section and mass distribution, and also found to be negligible (<0.01%). 
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FIG. 1; (Color online) Normalized combinatorial background (blue) for all and individual targets compared to data (black). 



Lepton pair production also has a background of random combinations of e+e~ pairs due to uncorrelated sources 
occurring within the same 2 ns CEBAF beam bucket. The most salient feature of the uncorrelated sources is that they 
produce the same-charge lepton pairs as well as oppositely charged pairs. The same-charge pairs (e+e^ and e~e~) 
provide a natural normalization of the uncorrelated background. This is also true for the measurement of opposite-sign 
pions or muons for which the combinatorial method has also been used in the past [sgI . [37| . This method has also 
been used in the extraction of resonance signals and proton femtoscopy of eA interactions [s^ . 

The combinatorial background is determined by an event-mixing technique. The electrons of a given event are 
combined with positrons of another event, as the two samples of electrons and positrons are completely uncorrelated. 
This produces the phase-space distribution where electrons and positrons are actually from different processes but 



5 



lying in the same event. 




e+e- Invariant Mass (GeV) 




FIG. 2: (Color online) Result of the fit to the e invariant mass obtained for the carbon (top) and Fe-Ti (bottom) data. 
The curves are calculations by the BUU model [40l. IZlj for various e^e~ channels. 

The mixed opposite-charged leptons chosen from samples of uncorrelated events were used to estimate the shape 
of the combinatorial background. This distribution was then normalized to the number of expected opposite-charge 
pairs. The result is shown in Fig. [T]for the individual targets. 

To simulate each physics process, a realistic model was employed and corrected for the CLAS acceptance. The 
events were generated using a code based on a semi-classical Boltzmann-Uehling-Uhlenbeck (BUU) transport model 
developed by the Giessen group that treats photon-nucleus reactions as a two-step process [42|. In the first step, 
the incoming photons react with a single nucleon, taking into account various nuclear effects, e.g. shadowing, Fermi 
motion, coUisional broadening, Pauli blocking, and nuclear binding. Then in the second step, the produced particles 
are propagated explicitly through the nucleus allowing for final-state interactions, governed by the semi-classical BUU 
transport equations. A rather complete treatment of the e~ pair production from 7A reactions at Jefferson Lab 
energies using this code can be found in Ref. 44]. 

The expected combinatorial background distributions are subtracted from the e+e^ effective mass distributions. 
The peaks of the lo and </) vector mesons are prominent in the invariant mass spectra, and one can determine the 
normalization of these narrow peaks rather easily. The shape of these peaks and the w Dalitz channel are well described 
by the BUU model where the ratio of a; to a; Dalitz decay was also fixed to their branching ratios. These distributions 
were fit to the data, then the resulting normalized heights were subtracted, leaving just the experimental spectra of 
the p mass. These fits for carbon and iron/titanium are shown in Fig. [2l 

The extracted p mass distributions are then fit with the exact spectral functional form obtained from calculating 
the cross section of production of the p meson including the Icptonic decay width [3, El] ■ The results of the 
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FIG. 3: The result of a simultaneous fit to the p mass spectra and the ratios for D2 (top), C (middle), and Fe-Ti (bottom). 




FIG. 4: Fits to the p mass distribution from D2 data using a Breit-Wigner function (dashed line) and Breit-Wigner//!^ (solid 
line) . 
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TABLE L The mass and width (in MeV) of the p meson obtained from the simultaneous fits to the mass spectra for each target 
and the ratio to Da compared to the result of the BUU simulations. The masses are consistent with the PDG values (770.0 ± 
0.8 MeV) and the widths are consistent with the coUisional broadening (150.7 ± 1.1 MeV). 



fits superimposed on the data are shown in Fig. O and the results are tabulated in Table [J The fits describe the 
data very well, and while the width of the p meson is consistent with the natural width of 150 MeV and collisional 
broadening [43|- that is also included in the BUU calculations- it is not compatible with the doubling of the p width 
reported by NA60 [HI. 

A good approximation for the p spectral function used to fit the data is a Breit-Wigner function divided by p^, 
where p is the mass of the e^e~ pair. The factor l/p^ comes from the photon propagator (l//i^) that couples to the 
p meson in the dilepton decay diagram, multiplied by a phase-space factor. Indeed the fits to the Breit-Wigner//z"^ 
rather than a simple Breit-Wigner function describe the data very well. For example, for the D2 target a simple 
Breit-Wigner fit gives a per degrees of freedom = 3.9 while a Breit-Wigner//j,'^ gives a psr degrees of freedom 
— 1.08 (see Fig. d]). The sensitivity of the fits to the 1/p^ factor indicates that the systematic uncertainties in the 
background subtraction are insignificant, and the p spectra are cleanly extracted. Similar results are obtained for the 
heavier targets, C and Fe-Ti, where the uncorrelated background is proportionally larger. 

Based on the notation of Ref. @ and Eq. [2 we obtain the shift parameter a — 0.02±0.02 for the Fe-Ti target with 
p momenta ranging from 0.8 to 3.0 GeV. This is consistent with no significant mass shift predicted by the calculations 
of Ref. 0, and those of Ref. [i^, at p vector meson momenta > 1 GeV. 

The total systematic uncertainty for the measured a due to various sources is estimated to be Aa = ±0.01 (3^ . 

Our result sets an upper limit of a = 0.04 with a 95% confidence level. This does not favor the prediction of 
Rcfs. 8] and 19] for a 20% mass shift and a — 0.16±0.06 respectively, and is significantly different from other similar 
experiments |2i, |2^, [sO], where a = 0.092±0.002, with no broadening in the width of the p meson. Our results are 
also not necessarily inconsistent with the result of the experiment in Ref. [3l| that measures a —14% shift in the mass 
of the uj meson, since different medium modification mechanisms are indeed expected for p and uj mesons (48l . l49j . 

The extracted experimental p mass spectrum with the unique characteristic of electromagnetic interactions in both 
the production and decay, is well described by the p functional form obtained from the exact calculations given 
in Refs. [3, \^ with no modification beyond standard nuclear many-body effects. With the availability of 
more sophisticated theoretical models and improved analysis techniques, future experiments with higher statistics are 
expected to make a conclusive statement about the momentum dependence of the in-medium modifications and the 
nature of the QCD vacuum. 
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de la Recherche Scientifique and Commissariat a I'Energie Atomique, the Korea Research Foundation, the U.K. 
Engineering and Physical Science Research Council. Jefferson Science Associates (JSA) operates the Thomas Jefferson 
National Accelerator Facility for the United States Department of Energy under contract DE-AC05-06OR23177. 
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